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ABSTRACT 


The purpose of this thesis is to generate a tactical decision aid (TDA) capable of 
calculating the probability of kill of a submarine when targeted with a vertical launched 
(VLA) anti-submarine rocket propelled torpedo (ASROC). In determining the submarine 
specific probability of kill (Pk), the Passive Contact Tracker and Kill Probability (PACT- 
AKP) TDA calculates the submarine’s position and its area of uncertainty (AOU) based 
on single or multiple ASW passive sensor bearing-to-target inputs. 

In determining the target’s position and AOU, PACT-AKP employs an extended 
Kalman Filter that uses MTST movement and measurement models. In calculating 
ASROC probability of kill, submarine specific torpedo specific effectiveness (TEFF) data 
collected from NUWC Newport was used to generate the Pk algorithm. 

We can conclude that PACT-AKP not only assists the ASW team with target 
motion analysis (TMA), but also provides the commander with a credible target 
probability of kill prior to the employment of VLA ASROC torpedoes as a deliberate 


attack weapon. 
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EXECUTIVE SUMMARY 


With limited weapon payload and increasing number of potential threat 
submarines, if the U.S. surface fleet is to contribute to maintaining maritime dominance 
the submarine threat must be met with the prudent and effective use of anti-submarine 
warfare (ASW) munitions. This thesis explores and solves a current and tactically critical 
ASW problem that plagues the surface navy fleet; the probability of kill of a submarine 
when engaged with a vertical launched (VLA) anti-submarine rocket propelled light 
weight torpedo (ASROC). This thesis develops the passive contact tracker and 
probability of kill (PACT-AKP), an EXCEL-based Microsoft Windows platform tactical 
decision aid (TDA) designed to aid in the passive tracking of a target and more 
importantly, to determine the target’s probability of kill when engaged with a VLA 
ASROC. 


PACT-AKP only requires the user to enter passive contact information normally 
recorded and processed during passive target motion analysis (TMA). Upon demand, 
PACT-AKP provides a graphical depiction of the submarine’s position, 2-sigma area of 
uncertainty (AOU) and ASROC probability of kill. In addition, PACT-AKP’s probability 
of kill feature is flexible enough so that it gives the user the ability to calculate the 


probability of kill without using the program’s passive contact tracker feature. 


Through use of Kalman filter theory and probability theory, this study provides a 
probability of kill for ASROC employment. Although numerous naval systems employ 
Kalman filters and its variations, what makes this thesis unique is the derivation of the 
probability of kill algorithm. Prior to this effort, there was no antecedent work tackling 
this problem. Anti-submarine warfare operations are complex and beset with uncertainty. 
Determining a submarine’s probability of kill when engaged with a VLA ASROC prior 
to the actual deliberate attack engagement provides the commander with the necessary 


knowledge needed to make better decisions concerning the expenditure of limited 


xiii 


shipboard ASW munitions. The probability of kill algorithm developed herein is unique 
and provides a foundation for the development of future Navy-wide ASW doctrine 


governing the use of ASROC munitions. 
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I. PACT-AKP 


A. INTRODUCTION 


With limited weapon payload and increasing number of submarines world wide, if 
the U.S. surface fleet is to continue to maintain maritime dominance the ASW threat must 


be met with the prudent and effective use of ASW munitions. 


This thesis explores and answers a current ASW problem facing the surface fleet, 
estimation of the probability of kill of a submarine when engaged with a vertical 


launched (VLA) anti-submarine rocket propelled light weight torpedo (ASROC). 


This thesis develops the passive contact tracker and probability of kill (PACT- 
AKP), a Microsoft Windows platform tactical decision aid (TDA) designed to aid in the 
passive tracking of a target as well as to determine the target’s probability of kill when 


engaged with a VLA ASROC. 
B. BACKGROUND 


When faced with an ASW threat, both the commander and the ASW team should 
be equipped with readily available decision aids capable of providing valuable and 
practical information. Passive tracking of a submarine using target motion analysis 
(TMA) is normally conducted by hand, however, this tedious method can be easily 
replaced or the process aided by a TDA. Moreover, the commander should also have the 
necessary information available to make an informed decision before engaging in the 


expenditure of the warship’s limited ASW munitions to neutralize the ASW threat. 


Because current inventories of VLA MK46 light weight torpedoes are low, 
determining the probability of kill for a hostile submarine is important. As defined by 
NTTP-3-21.33 Surface Torpedo Attack Tactics / Countermeasures / Evasion Manual, a 
deliberate attack is offensive in nature and it is a planned coordinated employment of 
firepower intended to destroy a hostile submarine. The most effective ASW weapon for 


the prevailing circumstance should be used. Due to the short employment ranges of 


1 


surface vessel torpedo tube (SVTT) launched torpedoes that requires maneuvering well 
within the lethality range of submarine torpedoes, VLAs should be used for deliberate 
attacks. Surface vessels equipped with MK-41 vertical launch systems and armed with 
ASROC torpedoes are the preferred surface platform for a deliberate attack due to their 


standoff weapon engagement ability. 


The surface fleet is challenged to realistically determine the probability of killing 
a submarine when engaged with a VLA ASROC. Despite popular belief, the “Pk” display 
on the OJ-452 console does not refer to probability of kill, but rather how well the 
underwater fire control system (UFCS) of the ANISQQ-89 integrated ASW system is 
tracking a target in terms of up-to-date bearing-to-target measurements. The fleet is 


currently without a method for assessing kill probability for the actual weapon. 


During a major fleet exercise Valiant Shield 2006 (VS06), simulated VLA 
employment greatly exceeded actual magazine capacity. It was observed that surface 
escorts followed a policy of “classify with ordinance” or using a weapon in the water to 
determine if the contact was an actual submarine (Naval Mine and Anti-submarine 
Warfare Command, (Forthcoming)). VS06 not only identified a weapon misuse issue but 
also a problem with kill probability determination regarding the VLA ASROC weapon 
system. The following year, exercise Valiant Shield 2007 (VS07) produced similar 
results. As published in TM 3-21.1-08 Vertical Launch ASROC (VLA) Employment, 
escorts expended more simulated torpedoes than physically possible. Based on 
observations from VS06 and VSO7, one of the logical conclusions is that commanders 
had to know how and when to use ASW weapons more prudently (NMAWC, 
(Forthcoming)). However, without the ability to ascertain the success of the weapon 
engagement as measured by the probability of kill, the commander must necessarily 


make less than optimal weapon expenditure decisions. 
C. PURPOSE 


This thesis’ purpose is to create a TDA capable of assisting the antisubmarine 


warfare team with the passive sensor tracking of an ASW threat and assist the 


commander in making prudent decisions regarding the employment of the VLA ASROC 
torpedoes by calculating the probability of kill (Pk). 


D. SCOPE AND METHODOLOGY 


PACT-AKP requires the user to manually enter passive sensor bearing-to-target 
data normally used while conducting TMA. PACT-AKP then processes these inputs 
using an Extended Kalman Filter (EKF) that employs a maneuvering target statistical 
tracker (MTST) movement and measurement model to estimate the position and velocity 
of the target. It graphically depicts the target’s area of uncertainty (AOU) as derived from 
the EKF process and calculates its probability of kill when engaged with either a MK-46 
or MK-54 ASROC. 


PACT-AKP uses unique submarine specific torpedo effectiveness (TEFF) data 
along with information derived from the EKF process to calculate the probability of kill. 
PACT-AKP’s kill probability algorithm uses TEFF data as determined by the Navy’s 
Weapon Analysis Facility (WAF) and published in TACD&E 00-20, 01-14, 02-11, 03-14 
and 06-XX . In this thesis, only unclassified data is used to demonstrate PACT-AKP. 


E. BENEFITS OF PACT-AKP 


Although the use of Kalman filters for the purpose of TMA is not a novelty, what 
sets PACT-AKP apart from any other TDA in its class is its fleet-wide accessibility, 
portability and above all, its unique feature of determining probability of kill. 


PACT-AKP enhances the capabilities of the Surface Navy by providing an 
effective passive ASW tracking process, and provides the commander with invaluable 
insight regarding VLA ASROC engagement allowing the commander to make more 


informed tactical decisions. 


In addition, PACT-AKP’s kill probability algorithm can easily be updated as new 
weapon effectiveness data is made available and its target tracking capability is scalable 
with slight modification of its embedded visual basic (VBA) code. PACT-AKP enhances 
the surface fleet’s ability to employ the VLA ASROC weapon system. 
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HW. PACT-AKP METHODOLOGY 


A. INTRODUCTION 


The purpose of this chapter is to provide the underlying theory and computational 
methodology used by PACT-AKP to both passively track and determine the kill 
probability of a target if engaged with a VLA ASROC torpedo. 


Throughout this chapter symbols for vectors and matrices will be bold, scalars 
will be italic and “<«—” is the replacement symbol. The superscript t on a matrix means 


“transpose.” 
B. METHODOLOGY 


PACT-AKP employs a variant of the Kalman filter to process user inputs in order 
to determine and project the state of the target from time ¢ to time ¢ + At. This variant is 


the extended Kalman filter (EKF). 


A basic Kalman filter is a method for recursively updating an estimate p of the 
state of a system by processing a succession of measurements Z. After each measurement 
step a new estimate is produced. In PACT-AKP, the measurement Z is bearing-to-target 
while p is an estimate of X, the location of the target in Cartesian coordinates. The 
associated covariance matrix X is used in calculating the dimensions of a 2-Sigma 
Ellipse area of uncertainty (AOU). The dimensions of the semi-major and semi-minor 
axis of the ellipse, as well as the estimated target position, along with TEFF data, are 


used to calculate kill probability. 
Cc. EXTENDED KALMAN FILTER 


In a basic Kalman filter, the measurement Z is related to X through the 
equationZ =HX+V. Z is the bearing-to-target measurement, H is the measurement 
matrix that describes how the bearing measurement depends on the target state and V is 


are independent Gaussian errors associated with the sensor. Because the relationship 


between the bearing measurement Z and the target state is non-linear, the matrix H is 
obtained by linearizing the non-linear function of Z The logic is as follows, if 
Z = f(X)+ V then the approximation Z = f(wWwt+HXK-p)+V, where 
H=df(X)/dX|X=n is the matrix of the first partial derivatives. The non-linear 
function of the state variables is approximated by the first order terms of a Taylor series 


expansion about p. Except for the fact that H depends on p, the approximation is a linear 


function of X (Washburn, 2007). 
1. Maneuvering Target Statistical Tracker (MTST) 


The Maneuvering Target Statistical Tracker was developed by Daniel H. Wagner 
Associates in the early 1980s and was chosen by the Navy as a Standard Tracker for at- 
sea targets. Currently, it is used in many tactical data systems, including the Global 
Command and Control System (GCCS) and Tomahawk Weapons Control System 
(TWCS). MTST builds upon the principles of Kalman filtering theory and makes the 
filtering process capable of processing non-linear measurements such as the ones 


associated with bearing-to-target measurements. 


In MTST, the state vector is a 4 x 1 matrix taking the following form: 


Long xX 
Lat ye 

= > 
Longvel Xvel 
Latvel Yve/ 


The extended Kalman filter used in PACT-AKP employs MTST movement and 
measurement models. Wagner (1989) explains in detail the mathematical theory behind 
each type of model. What makes the MTST models different from basic Kalman filtering 
are the particular formulations of the matrices, Q which are derived from the Ornstein- 
Uhlenbeck (O-U) process. The velocity components of the state vector are assumed to be 
independent Ornstein-Uhlenbeck processes (Wagner, 1989; Washburn, 2007). Both the 


movement model and measurement model are presented as applicable to PACT-AKP. 


a. PACT-AKP MTST Movement Model 


PACT-AKP’s MTST movement model has two embedded parameters that 
match submarine specific operating characteristics, and a user input parameter that 
represents the time in between measurements (Aft). The embedded parameters are the 
relaxation time (z) and the root mean squared (RMS) speed (s). These parameters are the 
fundamental building blocks of the MTST movement projection model and are embedded 
in PACT-AKP. 


The relaxation time (z) is the average time a target travels on a particular 
heading before changing course. PACT-AKP assumes 7 = 0.5 hrs. The RMS speed (s) is 
the average instantaneous speed of the submarine in knots. PACT-AKP assumes s = 7 
kts. 


To project the estimated state vector forward from time ¢ to time ¢ + At, 
MTST uses the updates p << ®p and X< O©L@'+Q. The matrices M and Q are given 
by: 


10 6, 0 c, 0 c, 0 

Of 1L20y 3B; [0 ¢, 0 ce, 
lo 0b 0] 16,0 «, 0 

000 b, 0 ¢, 0 ¢, 
Where: 


let c =.5s°, where s* has units of nm? /hrs” 
b, =(1—6,)r, where 6, has units of hrs. 





b, =exp'“”, where 0, is unitless. 
Cree (2Az —(3-4b, + b7 yr) tT, where ¢, has units of nm’. 
2 
C,=—, where c, has units of nm’ /hr. 
rE 


C,=C (1-b% ) where c, has units of nm’ /hr’. 


b. PACT-AKP MTST Measurement Model 


In the user interface of PACT-AKP, the user is required to enter bearing- 


to-target observations (Z) and the associated passive sensor error. Z is then modeled as: 
Z=HX+V 


In PACT-AKP, H is a 1 x 4 matrix and the V ’s are independent Gaussian 
vectors with mean zero and a | x | covariance matrix R. The matrix R is obtained by 
squaring the standard deviation of the measurement error of the bearing-to-target 


observation, as inputted by the user. The user is in all cases required to input the accuracy 


of each bearing, along with the bearing itself. 





Figure 1. Single Station Measure the Bearing to a Target Located at (X, Y).! 


Figure 1 provides a visual basis for the derivation of H. The location of 
the target is(X,Y)', with the observers (measuring stations / ships) located at known 


points (x, y). Let (D, @) be the true range and bearing from the observer to the submarine. 


Consider the measurement Z =9+V, where 9 = arctan Ged . The measurement is a 


X-x 


non-linear function of the state. Since: 


! All Figures and Tables within this thesis were created by the Author. 
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dO a cos(@) ae dO _ —sin(@) 


, and since @ does not depend on the 


dY D dX 








velocity component of X, H is [—sin(@)/D, cos(@)]/D, 0, 0] 


Because both @ and D depend on the unknown state of the target, H must 


be evaluated by inserting the latest estimates of @ and D. 


The target state is updated following a well defined sequence of 


calculations: 
1. Calculate the Kalman Gain: K = XH‘ | HEH' + R} 
2s Calculate p (post-observation state of the target): p<u+KS, where 


S =Z-8, and @ is the bearing to the expected target location. 


ae Calculate the post-observation covariance matrix: & <— [1 - KH] x, where I is the 
identity matrix. 


Wagner (1989) and Bailey (1992) further describe the EKF measurement 


model. 
Cc. PACT-AKP Initial Target Position Guess Entry 


PACT-AKP’s EKF requires a guess at the target’s state in order to get the 
filtering process started. The initial target position or “position guess” in PACT-AKP is 
entered by the user. The upper left 2 x 2 corner of the covariance of the target position is 
hard-coded into the program with the value 6.25 nm’, which indicates that PACT-AKP 
assumes the initial target position guess is an “educated guess” and accurate to within 
approximately 5nm (2Sqrt(6.25)= 5). This setting is PACT-AKP specific because the 
assumption is that when conducting passive TMA the operator will have a general idea of 
the bearing and distance of the target from the sensor by considering environmental 


conditions and sensor performance. 


Like the initialZ, the target position guess corresponds to p at time 0. 
Once the later is inputted or initialized, all subsequent calculations correspond to either 


the movements or measurements. 
d. PACT-AKP Coordinate System 


Most of the filtering process that takes place using MTST is performed on 
a flat earth. PACT-AKP requires the user to enter all inputs concerning measuring station 
position and target position guess in spherical coordinates of latitude and longitude, 


PACT-AKP later converts these coordinates into the X and Y coordinate system. 


The coordinate system used to describe the position of the measuring 
station as well as the target is the local East/North (Cartesian) coordinate system, which 
is defined at the point of the Earth that touches the tangent plane. The local East 


represents the X axis and the local North the Y axis. 


Almost all calculations are performed in the flat earth coordinate system 


with the exception of the bearing to the expected target location (0... ie), Which 1s 
calculated using spherical trigonometry. 0... ..q 18 the EKF’s best prediction of Z, based 


on all history previous to the 7” measurement (Washburn, 1982). PACT-AKP calculates 


O.. .cceqd USING spherical trigonometry in order to maintain measurement accuracy. For 


exp ecte 
instance, a one degree difference between two lines of bearing (LOB) to the same target, 
one LOB calculated in spherical coordinates and the other using flat earth coordinates can 


easily over a long distance generate a lateral difference of several miles. 
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e. MTST Shock and Dimensionless Shock 


exp ected 


S=Z-0 


expected 





Figure 2. Bearing-to-target Measurement vs. Expected Bearing-to-target. 


From Figure 2, the difference between what is actually measured (Z ) and 


O.xrecteq 18 called shock (S). From the measurement model Z = HX+V, the shock is the 
linear combination of the independent variables X and V. The approximate expected 
shock is 0 and the variance of S is HXH'+R, which is the denominator of the Kalman 
gain computation, and is the means for evaluating when the shock is too large. The 


dimensionless shock (DS) is given by: 


DS = S’(H=ZH‘ + RY" 


As Washburn (2007) explains, if the shock has m components, DS should 
be a scalar random variable that has a Chi-square distribution with m degrees of freedom. 
When the DS becomes large in comparison to m, it is usually because the filter has lost 


track. PACT-AKP displays DS after updating the state of the target. 


D. PROBABILITY OF KILL CALCULATIONS 


1. Target Relative to Weapon Water Entry Point AOU 


In calculating the target’s AOU, PACT-AKP uses the covariance matrix output 
from the EKF. The AOU calculated is a 2-Sigma ellipse which is an equiprobability 
contour that contains the target state with probability equal to 1- exp (-2) or .865 
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(Washburn, 1989). The ellipse has an inclination 2 and semi-major axis 2s; and a semi- 
minor axis 2s7 In calculating the dimensions of the ellipse we use the following matrix 
notation recalling that the output covariance matrix of the EKF is a 4 x 4 matrix 
consisting of 16 numbers of which we are only interested in using the upper 2 x 2 corner 


(four numbers). 





An ASROC will not hit precisely the spot at which it is aimed. The difference 
between the aim point and the water entry point is assumed to have a circular normal 
distribution with a circle error probable (ASROC CEP). In calculating the ASROC’s fire 
error (ASROC o), the ASROC CEP is divided by ./2In(2) (Washburn, 2002). PACT- 


AKP combines the ASROC o with both semi-major and semi-minor axis derived from 


the EKF covariance matrix as depicted below: 


S, = ,J(s? + ASROCo*), S, = ,|(s3 + ASROCo” ) 


The resulting semi-major and semi-minor axis are modified to account for the 
ASROC fire error, thus the target’s area of uncertainty becomes relative to the ASROC 
water entry point. The importance of this transformation becomes evident in the 
probability of kill calculation by allowing for the use of Washburn’s “EllipQ” VBA 
function. EllipQ builds upon Gilliland’s power series derivations for determining the 
probability of a hit within a specified target-centered circle when it is assumed that the 
causative missile guidance error is distributed according to the bivariate normal 


distribution (Gilliland, 1962). 
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Let R be the location of the target relative to the water entry point of the weapon. 
R is a critical random variable (RV) for the kill probability calculation and its cumulative 


density function (CDF) very difficult to calculate except for in the case where S, =S,,. In 


this case the CDF of R is: 





Pr(RSr)=1 exp 


Unfortunately, the problem at hand does not meet the above condition. The CDF 
of R as calculated by Washburn’s EllipQ function derived from Gilliland’s CDF of R 
derivation power series eq. 13. Gilliland’s method can be applied to end-point accuracy 
problems entailing conditions such as the one explored in this thesis: indefinite target 
location (Gilliland, 1962). Passive target tracking does not yield an exact target location, 
but rather a target position contained within a 2-Sigma Ellipse with 86.5% probability. 
Because of this uncertainty, the indefinite target location condition is applicable and the 
aforementioned method can be exploited with slight modification so that it is applicable 
to the problem at hand, the probability of kill of a submarine when engaged with a VLA 
ASROC. 


Washburn’s EllipQ function is a fundamental building block in the process of 
calculating the target’s probability of kill. EllipQ requires three inputs. The first input is 


the lethal radius of the weapon. The second and the third inputs are S, and S, as 


calculated for determining the target relative to weapon water entry point AOU. 
2. Torpedo Effectiveness (TEFF) 


Submarine type, evasion profile, and torpedo type and torpedo running profile 
specific TEFF data generated by the Weapon Analysis Facility (WAF) was provided by 
NUWC Newport. Because this information is classified SECRET, data similar in format 
was used to exercise the PACT-AKP’s Pk algorithm. A secret version of PACT-AKP that 
considers classified TEFF data can be made available. Figure 3 is a visual display of what 
TEFF data looks like prior to processing and serves as the anchor-point for further TEFF 


discussions. 
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Score Average Teff Zone (RangefAspect) 
60%. Splash-point= 7OOyds, Asp 113° 
600-S800yds, Aspect 90° to 136° 
30 400-1000¥ds, Aspect 68° to 158° 
Unsatisfactory =30% Everything Else 


Figure 3. TEFF Data Display. 














As observed from Figure 3, TEFF data is arranged according to angular (aspect) 
sectors defined by both their angular range and distance from center point. The center 
point (0, 0) of the graph represents the location of a particular submarine proceeding in 
the direction of the arrow. The sectors are divided into three categories: excellent, 
satisfactory, and unsatisfactory. Each sector is assigned a probability score (TEFF) that 
represents the drop zone torpedo probability of kill against the submarine. The 
probability of kill is based on which sector relative to the target the torpedo enters the 


water. 
3. Probability of Kill Calculation 


This thesis explores a unique, complex and tactically relevant problem that deals 
with target location uncertainty, weapon system fire error and a weapon system that 
cannot be described as having a lethal radius but rather a probability of kill as described 


in the TEFF section. 


Although the use of EKFs for tracking targets is not a novelty, PACT-AKP’s 
probability of kill algorithm is the only one of its kind. 
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a. Derivation and Assumptions of PACT-AKP Probability of Kill 
Algorithm 


From the TEFF data and as depicted in Figure 3, PACT-AKP uses the 


category sector distances (r) from center-point 7, =0 to the outer radius of the n” TEFF 


sector; n= 1,..,N. Where ris the outer radius of the mn” TEFF sector. 


n 


Let O, = EllipQ(r,, S,,S,). 


Let K be the event that the ASROC kills the submarine. Figure 3 shows 
P(K |r,@), where (7,0) is the location of target relative to the ASROC water entry point 
in polar coordinates. In spite of the recommended “splash point,” we assume that the 
ASROC is aimed directly at the submarine (the origin in Figure 3), and that the angle 0 is 


uniformly distributed [0, 27]. The logic behind these assumptions is as follows: 


(1) The submarine’s direction is unknown. Despite the best 
efforts of the TMA process, at the moment of weapon engagement it is plausible that the 
true direction of the submarine is unknown. The ASW team as well as the underwater fire 
control system may have a good idea of the general direction of the submarine, but one 
must not dismiss the possibility that the submarine may have changed course along the 
way since the last measurement. Because of this uncertainty, we assume that the 


likelihood of the submarine heading in any particular direction is the same. 


(2) The time lapsed between ASROC launch and water entry is 
minimal. Because of this assumption, leading the target will have little effect and it is 


assumed that the ASROC is much faster than the submarine. 


(3) As depicted in Figure 3, the area comprising the “excellent 
zone” only accounts for 5.7% of the total TEFF area depicted. The maximum value of the 


density of the location of the target relative to the ASROC water entry point is given 





(5 : } The probability of the ASROC hitting one of the excellent zones is the 
T 12 


product of the area of the excellent zone and the density. Assume that the target relative 
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to ASROC water entry point 2-Sigma ellipse AOU dimensions are S, = S, =500 yds , 
then the probability of hitting one of the excellent zones is only .07, even if one aims at 
the subject zone. 

Since @ is random, the kill probability depends only on r, the 
distance between the ASROC water entry point and the submarine. Let P represent the 


ring weighted TEFF score; n= 1,... N. 
PPK Ss [PK |r, 0)d8, for r,. <r <r, 
27 °0 


In Figure 3 for example, N= 4 and r, =400yds, r, = 600yds , 
r,=800yds and r,=1000yds. Also, let P,,, be the kill probability when the miss 
distance exceedsr,,. The kill probability in the third ring (n= 3 and 7, = 800 yards) and 


assuming that “< 30%” means 0.15 for the unsatisfactory region is: 


Pp EO) £5 Aeon eee a00yards 
16 4 





The above equation represents the weighted TEFF average for ring 


three defined by 600yrds <r <800yds. P,,, always corresponds to “unsatisfactory” in 


+1 


TEFF data, and we will always take it to be half of the given bound. Given Figure 3, PR 
would be .15. 


b. Probability of Kill Algorithm 


PACT-AKP’s probability of kill algorithm bases its computations on the 
theory of total probability (TTP), as well as the inputs from the weighted sector TEFF 
and Washburn’s EllipQ function. 


In following the previous notation, the Pk algorithm takes as 


input P, and Q,. 
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Pk= YPOe —Q,)+ Py aQn 


The difference (QO, ,—Q,) is the probability that the ASROC water entry 
point is in the n” ring. 

Although the algorithm is simple and elegant, the complexity resides in 
deriving the various elements that comprise the equation. 


Consider PACT-AKP probability of kill algorithm and the TEFF data in 
Figure 3. If one wished to calculate the Pk for the submarine represented in Figure 3, the 


following sequence of calculations would take place: 


P -(9) =.15, for 0<7r<400yds. 
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Pie ad =.225, for 400 <r < 600 yds. 
pees ae 158) = 3, for 600<r <800yds. 
p,-( 38)+.15@) 


ie = .225, for 800 <r <1000yds. 


P, =.15, where P, is the probability of kill beyond the 4” ring. 


Recall that EllipQ is a three parameter function that 


requires r, S, and S,as inputs. Where r is the radius of the applicable ring and 
both S, and S,are the semi-major and semi-minor axis of the target’s 2-Sigma ellipse 
AOU. Let S, = 1000 yds and S, = 500 yds. The Pk algorithm would take on the following 


form and values: 
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O, (400, 1000, 500) = .8548 


Q, (600, 1000, 500) = .7085 


Q, (800, 1000, 500) = .5526 


OQ, (1000, 1000, 500) = .4099 


Pk = FIL-Q)+P,Q,-9,)+ RQ; -Q,)+ P,Q, —9,) + BQ, 
Pk =.15(.1452) +.225(.1463) +.3(.1559) + .225(.1427) +.15(.4099) 


Pk =.1951 


All TEFF data is embedded, and S, and S, vary according to the EKF 


output. In addition, because all EKF calculations are performed in nautical miles, the 


input data where applicable is converted to nautical miles. 
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HW. PACT-AKP PERFORMANCE TESTING 


A. SCOPE 


The correct kill probability calculation is premised on the proper tracking of the 
target. As referenced in Chapter II, the change in the state of the target is not only 
depicted in terms of its position and velocity but also described by its 2-Sigma AOU. 
Because the dimensions of the target relative to weapon water entry point 2-Sigma AOU 
semi-major and semi-minor axis are direct inputs for Washburn’s EllipQ function, it is 


necessary to verify that such inputs are consistent with what the EKF is supposed to do. 


The purpose of this chapter is to verify the proper functioning of PACT-AKP’s 
EKF. A series of performance tests are conducted to test both the measurement and 


movement models as well as the proper graphing of the target’s 2-Sigma AOU. 


B. PACT-AKP MEASUREMENT MODEL PERFORMANCE TEST 


1. Target Position Guess Reliability Test 


The target’s initial position guess will remain constant at Cardinal Coordinates (0, 
0) as well as the area reference point at (35N/120E). As described in Chapter I, the 
initial hard-coded covariance matrix that describes the reliability of the “target position 
guess” is pre-established to be 5 nm. This means that if no measurements were taken and 
t remained at 0, the 2-sigma equiprobability contour describing the target state would be 


circular with a radius of S5nm. 


The aforementioned statement holds true and is displayed in Figure 4. 
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(hrs) — 0. 
- j Passive 
z Tracker 
-0-0 / 120-0-0 


1 Axis (Nm) = 5.0000 
Axi 


[AOU = 78.54 sq. nmi. 














Figure 4. | Target Position Guess AOU Display. 


2. Single Sensor Measurement Model Performance Test 


The purpose of this performance test is to exercise the measurement model of the 
EKF and verify the proper graphing of the 2-Sigma AOU. The target’s initial position 
guess remains constant at Cardinal Coordinates (0, 0) as well as the area reference point 


at (35N/120E). The sensor standard deviation (co) remains constant at 1.25 degrees. 


a. PACT-AKP Measurement Model Single Sensor and Single 
Bearing-to-target AOU 


(1) Case One: North-South Measurement 














Measurement Station Position Catlesian 
Coord 
Latitude Longitude x Y | Brng-to-tgt 0 Time 
Deg | Min | Sec | Deg | Min | Sec Deg Deg Hr 
35 10 00 120 00 00 0 10 180 123 0 



































Table 1. | North-South Single Sensor Bearing-to-target Measurement Station Data. 
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Based on the data from Table 1 and pre-established conditions for 


this performance test the following hypothesis is formulated: 


Case One Hypothesis: If PACT-AKP’s measurement model is 
processing the single bearing-to-target measurement correctly, the semi-major axis of the 
2-Sigma Ellipse AOU is larger along the bearing-to-target axis. With the sensor o being 
1.25 degrees or .02 radians and the distance to the target at 10 nm, one would expect the 


semi-minor axis to be approximately 10 *(2*.02) = .436 nm. 





i 00 
a a E- PACT-AKP TARGET AOU AND POSITION DISPLAY(nm) PK Passive 
z - ars Tracker 
Est. position = 35-0-0 / 1119-59-59 


}AQU Semi-major Axis (Nm) = 5.0000 














Figure 5. | North-South Single Sensor Bearing-to-target Target AOU. 


As observed from Figure 5 and consistent with data from Table 1, 
the measuring station is located at coordinate (0, 10) with a_bearing-to-target 
measurement of 180 degrees. The resulting 2-sigma AOU has the semi-major axis 
oriented along the bearing-to-target axis (Y-axis) with a dimension of 4.9 nm and a semi- 


minor axis of length .435. The hypothesis is not rejected. 


(2) Case Two: East-West Measurement. The same performance 
test and hypothesis as in case one applies to this case. This case uses measuring station 


data prescribed in Table 2. 
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Cartesian 


Measurement Station Position 
Coord 











Latitude Longitude Bmg-to-tgt | o Time 








Deg | Min | Sec | Deg | Min | Sec ESS Deg Deg Hr 














Gis) 00 00 119 47 45 270 125 0 



































Table 2. East-West Single Sensor Bearing-to-target Measurement Station Data. 


The location of the measuring station is at (10, 0) and the bearing- 


to-target measurement is 270 degrees. The results are depicted in Figure 6. 





5 59-59 / 1119-59-59 | 
|AOQU Semi-major Axis (Nm) = 5.0000 | 
AOU Semi-minor Axis (Nm) = 0.4347 
AOU = 6.83 sq. nmi. 














Figure 6. East-West Single Sensor Bearing-to-target Target AOU. 


From Figure 6, the semi-major axis is oriented along the bearing- 
to-target axis (X-axis). Because the same conditions regarding the sensor sigma as in case 
one have not changed, the semi-minor axis’ dimensions are proportionally equivalent. 


The hypothesis is not rejected. 
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b. PACT-AKP Measurement Model Single  Bearing-to-target 
Position Test 


Although the MTST measurement model performs all calculations on a 
flat earth it still accounts for the curvature of the earth when ascertaining the position of 
the target. If the measurement model is being implemented correctly this adjustment 
takes place at the Y axis coordinate of the target’s position. Table 3 depicts the measuring 


station data used to conduct this performance test. 




















Measurement Station Position ine 
Coordinates 
Latitude Longitude x Y Brng-to-tgt fo} Time 
Deg | Min | Sec | Deg | Min | Sec Deg Deg Hr 
35 00 00 119 00 00 | 49.1 0 270 5 0 



































Table 3. Target Position Test Measuring Station Data. 


After processing the data from Table 3, PACT-AKP estimates the position 
of the target to be located at (0, -.06). The negative coordinate in the Y axis is consistent 
with an adjustment for the earth’s curvature. Based on these results one can not reject the 


hypothesis that the measurement model is functioning correctly. 
3. Multiple Sensor Measurement Model Performance Test 


The purpose of this test is to verify that PACT-AKP’s measurement model is able 
to process more than one bearing-to-target measurement. The target’s initial position 
guess will remain constant at Cardinal Coordinates (0, 0) as well as the latitude / 
longitude reference point at (35N/120E). The sensor standard deviation (o) varies 


between 1.25 and 5 degrees. 


2D 





a. Two Sensor Measurement Model Test 


Table 4 depicts the measuring station data and Figure 7 the graphical 
display of the AOU. As it can be observed, the two measurements are orthogonal with 


the second measurement being less accurate than the first. 






































Station Measurement Station Position Canesten 
# Coord 
: : Pes o | Time 
Latitude Longitude tet 
Deg | Min | Sec | Deg | Min | Sec |X Y Deg Deg | Hr 
1 35 10 00 | 120 | 00 00 0 10 180 S) 0 
2 35 00 00 | 119 | 47 | 45 10 0 270 1,25 0 





























Table 4. | Multiple Sensors Measuring Station Data. 





PACT-AKP TARGET AOU AND POSITION a? Passive 
es Tracker 


AOU Semi-minor Axis (Nm) = 0.4347 
AOU = 2.25 sq. nmi. 














Figure 7. | Two Sensor Measurement Target AOU. 


From Figure 7, PACT-AKP estimates the position of the target to be at (0, 
-.01). In the case where no measurements were taken (Figure 4) the 2-sigma AOU had a 
radius of 5 nm, however, with two independent orthogonal measurements and different 


measurement accuracies the 2-sigma AOU becomes elliptical with a semi-major axis of 
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1.65 nm and semi-minor axis of .43 nm. The semi-major axis is oriented along the 
bearing-to-target measurement with the largest measurement error. On the other hand, if 
both bearing-to-target measurements shared the same sensor measurement error, the 2- 


sigma AOU would be circular with dimensions varying according to the sensor accuracy. 


C. PACT-AKP MOVEMENT MODEL PERFORMANCE TEST 


1. Change in Time Movement Model Test 


The purpose of this performance test is to ensure that the dimensions of the 
target’s AOU increase symmetrically with the passage of time given that no bearing-to- 
target measurements are taken. The target position guess is established at Cardinal 


Coordinates (0, 0) at start time 0. Table 5 depicts the results. 

















Time (Hr) Semi-major Axis (Nm) Semi-Minor Axis (Nm) 
0 5 5 
1 ieee 11.7 
2 157 15.7 
3 18.6 18.6 
4 pales os 














Table 5. | Change in Time 2-Sigma AOU Expansion. 
As observed from Table 5, the 2-sigma AOU expands symmetrically with time. 
D. PACT-AKP EKF PERFORMANCE TEST 


The purpose of this performance test is to exercise PACT-AKP’s EKF to ensure 
that it’s not only able to track a target, but also able to graph its corresponding AOU 
given multiple bearing-to-target measurements at different times with heterogeneous 


sensor bearing-to-target errors. This performance test is the most complex of all because 
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it involves simulating a moving target. For the purpose of this performance test, the 
generated target starts at (0, 0) and travel in a straight line from West to East at a speed of 
7kts. Actual Target track information is depicted in Table 6. In terms of measurement 
station bearing-to-target information, to ensure that PACT-AKP is able to track the target, 
three bearing-to-target cases are explored. In all three cases, PACT-AKP processes 
bearing-to-target measurements that modify the “true” bearing-to-target measurement by 
adding an error that is normally distributed with mean 0 and standard deviation of 5 


degrees. Measurement station data is depicted in Table 7. 
1. Composite EKF Performance Test 


True target location data is depicted in Table 6 and measurement station data is 


depicted in Table 7. 


























Target Position Cartesian 
ial Latitude Longitude Coomingtcs 
Hr Deg Min Sec Deg Min Sec x Y 
0 35 00 00 120 00 00 0 0 
1 35 00 00 119 51 20 - 0 
2 35 00 00 119 42 50 14 0 
3 35 00 00 119 34 15 21 0 
































Table 6. Actual Target Track Data. 


As depicted in Table 6, the target travels on a straight line at a speed of 7kts. At 
time 0, the target is located at origin (0, 0), but after a change of three time periods (3 hrs) 


it travels 21 nm eastward. 
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Cartesian 
Station # | Coordinates of the 
Measurement 
Station “Truth” | Case 1 | Case 2 | Case 3 o Time 

x ¥ Deg Deg Deg Deg Deg Hr 
1 0 10 180 180 185 182 5 0 
2 10 0 270 271 268 215 5 0 
1 7 10 180 178 176 182 5 1 
2 17 -10 He 316 311 315 5 1 
1 14 10 180 195 186 181 5 2 
2 24 -10 315 315 316 317 5 2 
1 21 10 180 180 181 172 5 1 
2 31 0 270 270 276 265 5 3 





























Table 7. | Measurement Station Data (Heterogeneous Sensor Error). 


Table 7 depicts the measuring station positions at the time bearing-to-target 
measurements are taken. Because there is no passive signal generated by the target the 
positions of the measuring stations were determined so that plausible bearing-to-target 
measurements could be generated given the known time dependent target position. For 
contrast purposes the column labeled “truth” represents the “true” bearing-to-target 


measurement. 


The end purpose of this test is to explore whether or not the various case scenario 
position estimate 2-sigma AOUs contain the true target. The true target position is 
depicted as a green diamond while the estimated target position is depicted as a smaller 


red diamond. 


2] 








Time Late (hrs) = 0.00 z 
Est. X = 0.00023317 PACT-AKP TARGET AOU AND POSITION DISPLAY(nm) Passive 
Est. Y = 0.14651732 Tracker 
Est. position = 35-0-8 / 119-59-59 
































AOU Semi-major Axis (Nm) = 1.6588 10 5 
AOU Semi-minor Axis (Nm) = 1.6372 
AOU = 8.53 sq. nmi. 
6 
4 
2 
eo 3 ae = fae me 2 ° = y 
2 
4 
6 
X-Lon (Nm) 
moe! 














Figure 8. | Case One Scenario at Time 0 After the First Two Measurements. 


Figure 8 displays case one target position estimate after processing the first two 
bearing-to-target measurements. In this case, the 2-sigma AOU contains the true target 


position. 
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Figure 9. | Case Two Scenario at Time 0 After the First Two Measurements. 
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Figure 9 depicts the 2-sigma AOU and estimated target position at time 0 for the 
case two scenario bearing-to-target measurements. In this case, the 2-sigma AOU 


contains the true target position. 










‘Time Late (hrs) = 0.00 
Est. X = -0.2828151 PACT-AKP TARGET AOU AND POSITION DISPLAY(nm) ~~ 
Est. 21559 


35-0-47 / 120-0-20 
|AO} Axis (Nm) = 1.7123 | 
|AQU Semi-minor Axis (Nm) = 1.6304 
AOU = 8.77 sq. nmi. 














Figure 10. | Case Three Scenario at Time 0 After the First Two Measurements. 


Figure 10 displays PACT-AKP’s target position and 2-sigma AOU estimate for 
case three scenario at time 0 after processing the first two bearing-to-target 


measurements. In this case, the 2-sigma AOU contains the true target position. 
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Figure 11. Case One Target Track at Time 1 After Four Bearing Measurements. 


Figure 11 displays the case one estimated target position after processing four 


measurements. In this case, the 2-sigma AOU contains the true target position. 
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Figure 12. Case Two Scenario at Time 1 After Four Bearing Measurements. 


Figure 12 depicts PACT-AKP’s target position and 2-sigma AOU estimate for the 
bearing-to-target measurement for case two after four measurements. In this case the 2- 


sigma AOU contains the true target position. 
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Figure 13. Case Three Scenario at Time 1 After Four Measurements. 


Figure 13 displays the target’s position and 2-sigma AOU estimation for the case 


three scenario at time 1. In this case the 2-sigma AOU contains the true target position. 
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Figure 14. Case One Scenario at Time 2 After Six Measurements. 


Figure 14 depicts PACT-AKP’s target position estimation for case one at time 


two. As depicted, the target’s true position is NOT contained by the 2-sigma AOU. 
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PACT-AKP processed a bearing-to-target measurement (195 degs) that we know is bad. 
At this point, PACT-AKP does not know any better since it has only processed six 
measurements. We can expect to notice a high dimensionless shock for one of the 


subsequent measurements. 








PACT-AKP TARGET AOU AND POSITION paca a 
. 1 Tracker 






Est. position = 35-1-12 / 119-44-12 
|AOU Semi-major Axis (Nm) = 3.5310 
AOU Semi-minor Axis (Nm) = 1.7906 
AOU = 19.86 sq. nmi. 

















Figure 15. Case Two Scenario at Time 2 After Six Measurements. 


Figure 15 depicts the case two scenario target position and 2-sigma AOU 
estimation after six bearing-to-target measurements. As observed, the 2-sigma AOU 


contains the target the true target’s position. 
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Figure 16. Case Three Scenario at Time 2 After Six Measurements. 


Figure 16 depicts PACT-AKP’s estimation of the target’s position and 2-sigma 
AOU when processing the case three scenario bearing-to-target measurements. In this 


case the 2-sigma AOU contains the true target position. 
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Figure 17. Case One Scenario at Time 3 After Six Measurements. 
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Figure 17 displays PACT-AKP’s estimation of the target’s position and 2-sigma 
AOU at time 3 for case one after six bearing measurements. The true target position is 


contained within the 2-sigma AOU. 
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Figure 18. Case Two Scenario at Time 3 After Six Measurements. 


Figure 18 depicts the target position and 2-sigma AOU estimation for case two 
scenario bearing-to-target measurements at time 3. In this case the 2-sigma AOU contains 


the true target position. 
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Figure 19. Case Three Scenario at Time 3 After Six Measurements. 


Figure 19 shows PACT-AKP’s target position and 2-sigma AOU estimation for 
case three scenario bearing-to-target measurements after six measurements. The 


generated 2-sigma AOU contains the target. 


In 11 out of 12 (91.6%) scenario cases, the generated 2-sigma AOU contains the 
“true” target. This phenomenon can be attributed to the fact that the initial position guess 
and the target’s initial position are the same. In addition, this test deals with a target that 
makes no course changes while PACT-AKP’s MTST movement model expects two 
course changes per hour. In theory, the true position should be in the AOU 86.5% of the 


time. 
2s Dimensionless Shock (DS) Test 


Another effective measure of EKF performance is evaluating the average 
dimensionless shock (DS) for all case scenarios. The average DS is expected to be one if 
the EKF is tracking a target that is behaving in accordance with the MTST movement 
model. Recall that PACT-AKP assumes the target makes two course changes in a one 
hour period. For this performance test, the simulated target is assumed to travel in a 
straight line (no course changes) from West to East for a period of three hours. According 
to PACT-AKP’s MTST movement model, PACT-AKP’s EKF expects the target to 
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change course no less than six times. By assuming the target travels on a straight line, we 
are not only lying to the EKF but surprising enough we are doing the EKF a favor. 
Because of the simulated target’s travel profile, one would expect the average DS to be 
less than one for all three scenarios. An average DS of less than one entails that PACT- 
AKP’s EKF is having no trouble tracking a target. The measurement dimensionless 


shocks and their averages for all case scenarios are depicted in Table 8. 





























Measurement Dimensionless Shock 
Measurement # Case One Case Two Case Three 
1 0 l 02 
2 0.004 .03 ell 
3 1.51 1.99 1.37 
4 0.08 01 .06 
5 0.7 1.27 2.37 
6 0.52 8 17 
7 3.311 2.14 2.65 
8 0.04 03 .06 
AVERAGE 0.77 8 85 

















Table 8. Measurement Dimensionless Shocks. 


Table 8 depicts the dimensionless shocks for all eight measurements processed in 
the sensor performance test. Dimensionless shock is an effective metric for ascertaining 
the accuracy of the bearing-to-target measurement in terms of what the EKF expects 
according to the conditions established in the MTST movement model. In essence, it’s a 
good indication of whether or not the EKF is maintaining track of the target. When 
processing a measurement the goal is for the DS to be as close to zero as possible. As 
observed, the average DS for scenario one is .77, which indicates that the EKF is doing a 
better job at tracking the target than the other two scenarios despite its measurement 
seven having the largest DS. Recall, that measurement five was a bad bearing-to-target 
input, and when the operator inputs the correct bearings for the next measuring period it 
temporarily throws the EKF a “curve ball.” PACT-AKP deemed measurement seven 
somewhat inaccurate as denoted by its DS. Scenarios two and three both have average 
dimensionless shocks of less than .9, which indicates that PACT-AKP’s EKF is on 


average having little difficulty tracking the target. 
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E. PACT-AKP PROBABILITY OF KILL ALGORITHM TEST 


The most unique feature of PACT-AKP is its kill probability function. As 
explained in Chapter II, the main inputs to the Pk algorithm from the EKF are the target 
relative to weapon water entry point AOU semi-major and semi-minor axis. Table 9 
depicts the various “target relative to weapon water entry point” AOU dimensions for the 
target state estimations calculated in the composite EKF case two scenario heterogeneous 


sensor error performance test. The weapon’s CEP is assumed to be 400 yards and the 


TEFF data used in this test is as shown in Figure 3 and annotated in Table 10. 





Target Relative to Weapon Water Entry Point AOU (Case Two Scenario) 

































































Measurement Time Semi-major Axis (Nm) Semi-minor Axis (Nm) 
0 1.81 1.61 
1 2.99 1.66 
Z 3.54 1.8 
3 1.9 1.49 
Table 9. | Case Two Scenario Target Relative to Weapon Water Entry Point AOU Data. 
ie Yards Nautical Miles Ve 
ra 400 975 - AUS 
t, 600 29624 PB, Dis) 
i 800 39499 ie a 
ve 1000 49374 2 225 
ie AS 
Table 10. | Converted TEFF Data. 


Table 10 shows the weighted average of the TEFF scores per region. The 


weighted TEFF score averages along with their distances from weapon water entry point 
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are the building blocks of the Pk algorithm. The target relative to ASROC water entry 
point 2-sigma AOU (heterogeneous sensor error) case two scenario data is depicted in 


Table 11. 




















Measurement Time | Semi-major Axis Semi-minor Axis AOU Area (Nm?) 
(t) “S$, (Nm) “5,” (Nm) 
0 1.81 1.61 9.19 
1 3 1.66 15.59 
2 3.54 1.8 19.97 
3 1.9 1.49 8.92 

















Table 11. Target relative to Torpedo Water Entry Point AOU Data. 


Based on the data from Table 10 and Table 11, the resulting measurement time 


specific Pk results are depicted in Table 12. 




















Measurement Time (t) Probability of Kill 
0 1534 
1 152 
2 1516 
3 51535 











Table 12. Pk Results. 


The probability of kill results as depicted in Table 12, show that the Pk against the 
submarine is at its highest at time 3 when the 2-sigma AOU area is at its smallest. These 
results show a relationship between the area of the 2-sigma AOU and Pk value 
highlighting the importance of not only the sensor accuracy, but also the frequency and 
number of bearing measurements. 
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In addition, suppose that S, = S, for values ranging from S, = S,=0to S, = S,= 
.5 nm. Recall that the CEP of the weapon is 400 yards or .197 nm. The Pk results are 
depicted in Table 13 and Figure 20. 

































































S, = S, (Nm) Probability of Kill 

0 .1976 
025 1983 
.05 .2003 
075 2031 
AI 2059 
h2d .2083 
15 .2097 
lds 2101 
.197 (400 yards) 2095 
22o 2077 
as 2054 
2h 2028 
Se, 1998 
a2 .1968 
oD .1939 
Bt 191 
4 1882 
425 1856 
45 1832 
475 181 
BS) 1789 











Table 13. Circular 2-Sigma AOU Probability of Kill Results. 


Table 13 shows that the submarine’s probability of kill is at its highest when S, = 


S,=.175 or 355 yards. The same holds true when S, = S,= .165 or 334 yards. 
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Figure 20. Pk vs. SJ. 


As observed from Figure 20, a small 2-sigma AOU may be in some cases 
detrimental because the ASROC will land too close to the origin, likewise, if the area of 
uncertainty is too large the value of the Pk will be low. In figure 20, a maximum Pk of 
.2101 is achieved when the 2-sigma area of uncertainty is between 173 and 195 yards?. 
As grim as this result may be, one must keep in mind that torpedo TEFF data varies 
according to the type and evasion profile of the submarine as well as the type of torpedo 
and torpedo run settings. Depending on the torpedo effectiveness scores associated with 


the various rings the Pk will sometimes be higher and other times lower. 
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IV. CONCLUSION 


A. INTRODUCTION 


This thesis explores and answers a critical tactical decision aid requirement that is 
currently required by the surface Navy fleet; the probability of kill for VLA ASROC 
torpedo launch. PACT-AKP makes use of Kalman filter theory. Numerous other naval 
systems employ Kalman filters and variations of such, but what makes this thesis unique 
is the derivation of the probability of kill algorithm. Anti-submarine warfare operations 
are complex and beset with uncertainty. ASW operators have to not only contend with a 
dynamic battle space, but also with a myriad of other factors that often make anti- 
submarine warfare operations overwhelming exercises of educated guessing. The work 
developed in this thesis aids in bounding this uncertainty by tackling and solving the 
uncertainty surrounding the probability of kill for VLA ASROC torpedo launch. The 
passive contact tracker and probability of kill TDA developed in this thesis provides the 
commander with a tool capable of assessing the ASW deliberate attack engagement 
success, thus, providing him the knowledge necessary to better make decisions 


concerning limited shipboard ASW munitions expenditure. 
B. RECOMMENDATIONS FOR FUTURE AREAS OF STUDY 
PACT-AKP does not address optimal weapon load-out problems. The following 


ASW relevant area of future study is recommended: 


Ls Implement PACT-AKP calculated Pk to determine combatant commander 
area of responsibility (AOR) VLA ASROC inventory requirements through 


campaign modeling. 


2 Explore optimal shipboard VLA ASROC inventory using Pk values and 


expected mission assignment. 
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3. Although PACT-AKP bases its computations of Pk on the assumption that the 
ASROC is aimed directly at the target, improvements may be possible by 


choosing other aimpoints. 
C. CONCLUSIONS 


PACT-AKP is an innovative and unique tactical decision aid that helps the 
commander make better informed decisions regarding the expenditure of limited ASW 
munitions. PACT-AKP enhances the TMA process by providing a mathematically sound 
platform from which the ASW Officer and commander alike can individually track a 
submarine based on passive bearing-to-target measurements, and compare the tracking 


results to the pencil and paper target motion analysis performed by the ASW team. 


Overall, PACT-AKP not only addresses a tactically ASW relevant problem, but 
also in the process of doing so, it provides the commander with an alternative method of 
performing TMA and target kill probability knowledge before the expenditure of VLA 
ASROC munitions. The probability of kill algorithm developed herein is unique and will 
provide a credible foundation for the development of future Navy-wide ASW doctrine 


governing the use of ASROC munitions. 
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